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The meson-cloud model of the nucleon consisting of a system of three valence 
quarks surrounded by a meson cloud is applied to study the electroweak structure 
of the proton and neutron. Light-cone wavefunctions are derived for the dressed 
nucleon as pictured to be part of the time a bare nucleon and part of the time a 
baryon-meson system. Configurations are considered where the baryon can be a 
nucleon or a A and the meson can be a pion as well as a vector meson such as the 
p or the uj. An overall good description of the electroweak form factors is obtained. 
The contribution of the meson cloud is small and only significant at low Q^. Mixed- 
symmetry S"-wave components in the wavefunction are most important to reproduce 
the neutron electric form factor. Charge and magnetization densities are deduced 
as a function of both the radial distance from the nucleon center and the transverse 
distance with respect to the direction of the three-momentum transfer. In the latter 
case a central negative charge is found for the neutron. The up and down quark 
distributions associated with the Fourier transform of the axial form factor have 
opposite sign, with the consequence that the probability to find an up (down) quark 
with positive helicity is maximal when it is (anti)aligned with the proton helicity. 

PACS numbers: 12.39.-x, 13.40. Gp, 14.20.Dh 



I. INTRODUCTION 



Since the discovery of the proton finite size by Hofstadter and coworkers more than 50 
years ago {l| , electromagnetic form factors have played aprivileged role in the investigation 
of the nucleon structure (for recent reviews, see 2|, [sl, Ij, [sl, ISj). 

In principle, these quantities reflect the strong interaction between quarks inside the 



2 



nucleon and should be described by quantum chromodynamics (QCD). In recent years im- 



portant progress has been made with lattice QCD simulations la, y, y, Q, Generic 
features of the baryon octet mass spectrum are reproduced well in quenched lattice QCD 
simulations (in which sea-quark contributions are neglected) and electromagnetic properties 
can be studied with pion masses as low as 0.3 GeV Qualitative agreement with the 



experimental data has been obtained 



e.g. the flavour dependence of the Dirac form 



factor Fi(Q^) 3|. The isovector nucleon form factors were calculated in both the quenched 
and unquenched approximations with configurations for pion masses down to 380 MeV j^. 
Small unquenching effects and results larger than the experimental data were found. Eval- 
uating the axial form factors as well as the irNN and vrA^A form factors it was possible 

n 

to check the Goldberger-Treiman relations p^. Qualitatively consistent results with ex- 
periment were obtained for the isovector form factor ratio Gp{Q'^)/Ga{Q'^) of the nucleon 
axial vector form factors At present, the state of the art is still limited by systematic 
errors which are related to the fact that calculations are performed on finite volumes, at 
finite lattice spacings and at quark masses which are still relatively large. In addition, the 
extrapolation to the chiral limit with the help of the chiral effective field theory requires 



15|, 



IGj (see also [171]) in order to 



calculations of higher order than presently available p^, 
account for the non-analytic behaviour of the form factors on the quark masses 181 ]. 

Therefore, model calculations are still a valuable tool for understanding the internal 
dynamics of the nucleon. In particular, many theoretical calculations have been done to 
investigate possible interpretations of the decreasing ratio of the proton electric to magnetic 
brm factors, UpG^^iQ"^) / G^iiQ"^) ^ as a function of (see, e.g., [l9| for an overview and 
20, l21 \ for the two-photon physics in elastic electron scattering) and we refer the reader to 
Refs. % Isj for a discussion of the different interpretations of nucleon electromagnetic form 
factors proposed in the literature. 

Here we want to consider the possibility that the physical nucleon is a bare nucleon 
surrounded by a meson cloud as a consequence of the spontaneously-broken chiral symmetry. 
As first discussed in the context of deep inelastic scattering 22, 1231] , a pion cloud can give an 
explanation of the flavor-symmetry violation in the sea-quark distributions of the nucleon 
thus accounting for the violation of the Gottfried sum rule [2J]. This cloud will manifest 
itself as an extension of the charge distribution of protons and neutrons, which should be 
observable in the electromagnetic form factors at relatively small values of . In fact, the 
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neutron charge density extracted from available data shows a positive core surrounded by a 
negative surface charge, peaking at just below 1 fm, which can be attributed to a negative 
pion cloud 25||. This is confirmed by the analysis of Ref. 26(] showing a pronounced bump 
structure in the neutron electric charge form factor G^{Q'^) (and a dip in the other nucleon 
form factors) around = 0.2 — 0.3 GeV^, which can be interpreted as a signature of a very 
long-range contribution of the pion cloud extending out to 2 fm. It must be said, however. 



that from dispersion relation analysis 27|, 1281] the pion cloud should peak much more inside 



the nucleon, at ~ 0.3 fm, and the desired bump-dip structure of Ref. [26| can only be achieved 
at the cost of low-mass poles close to the u mass in the isoscalar channel and to the three- 
pion threshold in the isovector channel 29|]. In addition, while confirming the long-range 
positively (negatively) charged component of the proton (neutron) charge density, a recent 
model-independent analysis of the infinite-momentum-frame charge density of partons in 
the transverse plane 3^ is suggesting that the neutron parton charge density is negative at 
the center. 

Mesonic degrees of freedom are naturally taken into account in the baryon chiral per- 
turbation theory that is the effective field theory of the standard model at low energies 
and small momentum transfer. The electromagnetic form factors of the nucleon have been 
calculated to fourth order (one-loop) in baryon chiral perturbation theory within the mani- 
festly Lorentz-invariant infrared regularization approach 3l|] and the extended on-mass-shell 



renormalization scheme 



32 



33| . The inclusion of vector mesons as explicit degrees of free- 
dom results in a considerably improved description, accurate up to — 0.4 GeV^. 



The problem of considering the meson cloud surrounding a system o 



has been addressed already in the past in a variety of papers (see, e.g., [34 



three valence quarks 
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36 



40 



41 



42| and references therein). Along the lines originally proposed in Refs. 
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3t 

44, 



39, 
45| 



and developed in [35|, in this paper a baryon-meson Fock-state expansion is used to construct 
the state \N) of the physical nucleon. In the one-meson approximation the state |A^) is 
pictured as being part of the time a bare nucleon, |A^), and part of the time a baryon-meson 
system, \BM). The bare nucleon is formed by three valence quarks identified as constituent 
quarks according to the ideas discussed in 46| . The model was revisited in Ref. 47|] to study 
generalized parton distributions where the meson cloud gives an essential contribution in the 
socalled ERBL region. We apply here the model to calculate the electroweak form factors 
of the nucleon. The baryon-meson system is assumed to include configurations where the 
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baryon can be a nucleoli or a Delta and the meson can be a pion as well as a vector meson 
such as the p or the u. 

In Sect, inithe relevant formulae of the meson-cloud model are collected and the light-cone 
wavefunctions derived. The calculation of the electroweak form factors of the proton and 
neutron is illustrated in Sect. IIIIl and the results are presented and discussed in Sect. HVl 
Concluding remarks are collected in the final section. Technical details necessary to calculate 
the vertex functions describing the transition to a baryon-meson state with vector mesons 
are given in Appendix S 

II. LIGHT-CONE WAVEFUNCTION OF THE NUCLEON IN THE MESON 

CLOUD MODEL 

The derivation of the nucleon light-cone wavefunction (LCWF) in the meson cloud model 
has been already discussed in Ref. ^j^- In this section we review some pertinent formulae 
necessary for the calculation of the form factors. In the meson cloud model the nucleon 
is viewed as a quark core, termed the bare nucleon, surrounded by a meson cloud. The 
mesonic effects are treated perturbatively, by truncating the Fock-space expansion of the 
nucleon state to the dominant components given by the bare nucleon and the state containing 
a virtual meson with a recoiling baryon. The corresponding quantum state of the physical 
nucleon (A^), with four-momentum = (padPtV) Pa^-l) = {VniVn) and hehcity A, can be 
written as 

Ip., A; N) = Vz\p,, A; AT) + E / E <P'A'''''\y^ k,) 

X \yp%, + yp^^, A'; B) \{1 - y)p+, -k^ + (1 - y)p^^. A"; M), (1) 

where the function (j)y^y,'^^'^\y,'k±) is the probability amplitude to find a physical nucleon 
with helicity A in a state consisting of a virtual baryon B and a virtual meson M, with 
the baryon having helicity A', longitudinal momentum fraction y and transverse momentum 
kx, and the meson having helicity A", longitudinal momentum fraction 1 — y and transverse 
momentum — k^. From the normalization condition of the nucleon state 

iP'^, A'; N\pN, A; N) = 2(27r)3p+5(p'+ - pt,)S^'\p'^^ - Pn±)Sxx', (2) 
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one obtains the following condition on the renormalization factor Z 



B,M 



with 



BM/N 



BM/N 7 



.l/2(N,BM) 
^A'A" 



(3) 



(4) 



From the definitions in Eqs. ([3]) and (HI), one can interpret the factor Z as the probability of 
finding a bare nucleon in the physical nucleon, and Pbm/n as the probability of fluctuation 
of the nucleon in a baryon-meson state. 

The probability amplitude (pyy''^^^ can be calculated using time-ordered perturbation 
theory in the infinite- momentum frame as explained in Ref. 47|. The final result reads 

1 



V,\AN,BM) 



(5) 



/y(l - y) Mf, 

where Vy y,{N , BM) is the vertex function describing the transition of the nucleon into a 
baryon-meson state, with squared invariant mass 

Ml + k2 Ml + k2 



y 



The vertex function V^, y, {N, EM) has the following general expression 48] 



Vly,{N, BM) = UNaipN, \K^^Xp{pM, X")^,{pB, A' 



(6) 



(7) 



where un- is the nucleon spinor, x and ip are the field operators of the intermediate meson 
and baryon, respectively, and a, (3, 7 are bi-spinor and/or vector indices depending on the 
representation used for particles of given type. The explicit expressions for the vrA^ and 
ttA cases have been derived in Appendix C of Ref. 47|, while the corresponding results 
for transitions with vector mesons are worked out in Appendix |X1 Because of the extended 
structure of the hadrons involved, one has also to multiply the coupling constant for pointlike 
particles in the interaction operator by phenomenological vertex form factors. These 
form factors parametrize the unknown microscopic effects at the vertex and have to obey 
the constraint F^BMiv, kj_) = Fnbm{^ — ?/, A:^ ) to ensure basic properties like charge and 
momentum conservation simultaneously 
form 

FNBA4{y,kj_) 



491] . To this aim we will use the following functional 

(8) 



exp 



/1#2 _ T\yT2 



2A 



BM 
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where A^^v/ is a cut-off parameter. 

For the hadron states of the bare nucleon and baryon-meson components in Eq. ([1]) we 
adopt a hght-cone constituent quark model, by using the minimal Fock-state wavefunction 
in the light-cone formalism, i.e. 



\ph,X;H)=Y: J 



1 


dx 








N 



[d'k^]iv^A 



N 



({xi, k_Li; Ai, ri}i=i,„,,Ar) H \^iPH^ P*-L) K n), (9) 

1=1 



where '^x {{xi,\<i±i]\i^Ti}i=i^,,^N) is the momentum LCWF which gives the probability 
amplitude for finding in the hadron partons with momenta (xjp^,pj_L = kj_L -|- Xiipni), 
and spin and isospin variables Aj and Xj, respectively. In Eq. ([9]) and in the following formulae, 
the integration measures are defined by 

dx 

. TV 



[d^kj 



\N 



n 



na)'('-H 

/ N 

2(27r)3d^kj 



(10) 
(11) 



v.=i2(2vr)^ , 
where the number of valence partons is iV = 3 and N = 2 for the baryon and meson 

case, respectively. As explained in Ref. 50|, the wavefunction ^P^''-^^ can be obtained by 

transforming the ordinary equal-time (instant-form) wavefunction in the rest frame into that 

in the light-front dynamics, by taking into account relativistic effects such as the Melosh- 

Wigner rotation, i.e. 



''■''^ ({^i' k_Li; Aj, Ti}i-i^„^N) 



N 



2(27r) 



'Mq fj-^ \Xi 



N 



X H ({k h=i,...,iv)n^;:i*(^c/(A:.)), 

^l,...,/ijV 



i=l 



(12) 



where vE^^''*^' is the canonical wavefunction, and D^J^^*{R(.f(ki)) are the Melosh rotations 
defined in Ref. |5o|. In Eq. (IT^ . uoi = -^m? + kf is the energy of the i-th quark, and 
^0 = Hi is the free mass of the system of A^ non-interacting quarks. 

In our model calculation, we take into account the meson cloud contribution correspond- 
ing to vr, p, and cj, with the accompanying baryon in the \BM) component of the dressed 
nucleon being a nucleon or a A. The instant-form wavefunction is constructed as the product 
of a momentum wavefunction, which is spherically symmetric and invariant under permu- 
tations, and a spin-isospin wavefunction, which is uniquely determined by SU(6)-symmetry 
requirements. 
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5l||, which 



In the case of the nucleon, we adopt the momentum wavefunction of Ref. 
reads 

- JJ^iT^, (13) 

with A^' a normahzation factor. In Eq. f[T^ . the scale (3, the parameter 7 for the power-law 
behaviour, and the quark mass m are taken as free parameters, and will be determined by 
a comparison with experimental data as explained in Sect. IIV[ 

The A is described as a state of isospin T = 3/2 obtained as a pure spin- flip excitation of 
the nucleon, with the corresponding momentum wavefunction equal to that of the nucleon 
in Eq. 



Furthermore, the canonical wavefunction of the pion is taken from Ref. 52| and reads 

r'^'KhM) = :;^y^,^^P[-kV{2a% (14) 

with k = ki = —k2, and the two parameters a = 0.3659 GeV and = 0.22 GeV fitted to 
the pion form factor data. The pha se of the pion wavefunction (IT^ is consistent with that 
of the antiquark spinors of Ref. 53 1 . 

The wavefunction of the p differs from the pion only in the spin component, with the 
rest-frame spin states of the qq pair coupled to J = 1 instead of J = 0. Similarly, the uj 
is described by the same spin and momentum wavefunction as the p, but with the isospin 
component corresponding to a singlet state. This choice corresponds to assuming an ideal 
mixing in the vector sector, since the effects of the (p — u mixing are irrelevant in the 
calculation of the meson cloud contribution to the nucleon form factors. For the same 
reason, also the effects of the — uo mixing are neglected. 

Finally, we need to specify the parameters entering in the vertex functions. The cutoff 
Kbm in Eq. ([8]) should in principle be different for each BM component. However, the Jiilich 



group 



541 ] and Zoller [4J] used high-energy particle production data to determine all the Kbm 



of interest, and found that the data could be described by two parameters: Ai for octet 
baryons and pseudoscalar and vector mesons, and A2 for decuplet baryons. We have chosen 
the values Ai = 0.61 GeV and A2 = 0.81 GeV, which are consistent with the ones adopted 
in the cloudy bag model 49|] to obtain a good fit to both the violation of the Gottfried sum 
rule and the measured sea quark contribution in the unpolarized parton distribution. For 
the NBM coupling constants at the interaction vertex we used the numerical values given 
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Table I: Coupling constants at the NBM interaction vertex. 



SnNtt 
in 


JNAtv 


2 

3nnp 


Innp 


J NAp 
An 


An^ 


13.6 


11.08 GeV-2 


0.84 




20.45 GeV-2 


8.1 



in Refs. |55|, l56j in the case of the n and the p. Instead, for the NNu coupling, Qnnlj, we 
used the result from the analysis of Ref. jsT] about the uj contribution to the unpolarized 
anti-quark distributions, which favors a much smaller value for Qnnu) than the one used to 
describe the nucleon scattering data. The numerical values for each of the BM states are 
summarized in Table [B 

Finally, with the specified parameters, the probabilities for each of the BM component 
in the dressed nucleon are 



Pntt/p 


Ppn'^ Ip 


~l~ Pnn+ Ip 


'^Ppn^'' Ip 


5.1%, 




Pa-k/p 


= PA++n 


-/p + Pa+ 


n"lp + PaOtt 


f/p = 2PA++7r-/p 


= 3.40%, 


Pnp/p 


= Ppp^'/p 


~l~ Pnp'^/p 


— '^Ppp^lp = 


0.11%, 




Pap/p 


= Pa++p 


-Ip + 


pO/p + PaOp+ 


Ip = 2Pa++p-/p = 


= 0.67%, 


Pnuj/p 


Ppui/p 


= 0.013% 









III. ELECTROWEAK FORM FACTORS OF THE NUCLEON IN THE 

MESON-CLOUD MODEL 



The Dirac and Pauli form factors Fi{Q'^) and F2{Q'^) of the nucleon are given by the spin 

'+ 

V 

= {P + q,l\4\p,l), (15) 

(16) 

^ ^ 

where = —q^. As was first shown by Drell and Yan 58], the calculation of the form 
factors is conveniently done in a coordinate frame where = 0. In particular, we will use 
a symmetric frame where the nucleon momenta are given by 

^Ir + qi/4 + 1 _ rM2, + qi/4 _ 



conserving and the spin-flip matrix elements of the vector current Jy = Jy + Jy 

{p + Q,^\Jv\P,\), 
{q. + iqy)F2{Q^) = 2M^(:p + g,-^|J+|p,^), 



Pn 



Pn 



Pn 



q± 

' 2 



Pn 



,Pn 
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p'n 



Pn 



(17) 



With such a choice, the processes with vacuum pair production are suppressed, and the 
current matrix elements can be computed as a simple overlap of Fock-space wavefunctions, 
with all off-diagonal terms involving pair production or annihilation by the current or vacuum 
vanishing. In the present meson-cloud model, we need to consider the contributions from 
the diagonal overlap between the bare-nucleon state, on one side, and the BM components, 
on the other side. Furthermore, the electromagnetic current is a sum of one-body currents, 
J+ = Y.B,M Jb + J Ml which involves individual hadrons one at a time. This corresponds 
to assuming that there are no interactions among the particles in a multiparticle Fock state 
during the interaction with the photon. Therefore the external probe can scatter either on 
the bare nucleon, |A^), or one of the constituents of the higher Fock states, \BM). As a 
result, the matrix elements of the electromagnetic current can be written as the sum of the 
following two contributions 



(18) 



In Eq. (1181) . is the contribution from the bare nucleon corresponding to the diagram (a) 
in Fig. [H and 61 is the contribution from the BM Fock components of the physical nucleon. 
This last term can further be split into two contributions, with the active particle being the 
baryon or the meson {SI(.m'm)b^^ ^ ^ 



B,B'M M'M,B 



{M'M)B 

jV'AiV 



(19) 




Figure 1: Electromagnetic interaction vertex for a bare nucleon (a), and virtual baryon (b) and 
meson (c) components of a dressed nucleon. 
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The 51^^ term in Eq. f|T9l) is schematically represented in Fig. [^b) and is explicitly 
given by 

B,B',M X,\',X" ^\^'^ ) 

where ks_L = Pb± - (1 - l/B)q±/2 and k'B'± = Pb± + (1 - 2/B)q±/2. 

Analogously, the contribution from the meson in the BM fluctuation is described by the 
diagram (c) in Fig. 1 and reads 

B,M,M' X,X',X" •' ^l^^J 

^^X^^iN.BM)^^ - VM, -kMx) - 1/M, -k'Ma)]*, (21) 

with = Pa/± - (1 - l/j\/)q±/2 and k'Af'± = Pa/± + (1 - yAf)q±/2. 

As a result, the contribution from the BM components in Eqs. fl2Ul) and fl^Tl) are obtained 
by folding the current matrix elements of the baryon and meson constituents with the 
probability amphtudes describing the distributions of these constituents in the dressed initial 
and flnal nucleon. In general, the current matrix elements (p + q\J^\p) appearing inside 
the integrals depend on the internal momentum of the baryon-meson state. However, as 
discussed in Ref. jsol, the kinematical nature of a light-front boost allows us to transform 
these matrix elements to a frame with p = 0, with the result 

h'x{Q^) = {p + qA'\J^\PA) = (M,q^,AV+|M,Ox,A). (22) 

As a consequence, the current matrix elements in Eqs. (1201) and (12 ip factor out of the internal 
momentum integration, and one flnds 

E E ix^x'm 

B,B',M A, A', A" 

X / dysj ||^4^:,'"'»"(,,.k«j iA;"''""'(,«.ki. (23) 

B,M,M' X,X',X" 

X. [a /■ d^PA/± ,Ajv(Ar,BA./)/-. . . r X'j^ {N,BM') , , .-.^ 

X J ^Vm J 2(27r)3 '^^"^ ^ ~ t*^^"^' " y^'' -Ki'±)] ■ 

(24) 



^AB'B)M 



(M'M)B 
Ajv 
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We also note that the sum in Eq. fll9l) over all the possible BM configurations leads to 
contributions from both the diagonal current matrix elements with the same hadrons in the 
initial and final state {B' = B and M' = M in Eq. (!23|) and (p^ . respectively), and the 
current matrix elements involving the electromagnetic transition between different hadron 
states (i.e. the terms with B' ^ B and M' M in Eq. ([23]) and ([24]), respectively). 



Finally, the current matrix element for the bare 
integrals of the hadron LCWF (see, e.g., Ref. [53, 



ladron states can be calculated as overlap 



6Q|) 



{p'h,,X'h',H'\4\ph,Xh,H) 

(25) 



A, 



where k'j_j = k±j + {1 — Xj)q± for the struck quark, and k_Lj = Xiq± for the spectator quarks. 

The convolution formulae derived for the electromagnetic form factors can be extended 
to the calculation of the proton axial form factor G'a(Q^)- In this case we need to consider 
the diagonal matrix element of the axial current J^, i.e. 

= + (26) 

The structure of Eqs. (ITS]) . (fT9]) . (!23l) and (12^ apply also for the axial matrix elements, 
with the difference that the active mesons in the 5liM'M)B contribution can only be vec- 
tor mesons. Furthermore, the LCWF overlap representation of the axial matrix elements 
between bare hadron states reads 

{p'h,A'h'.H'\J\\PhAh.H) 

= J2^jI1 sign(A,) / H[d2k^][vpJ' [^l({x„ k'^,; A„ T,},=,i,..,N)r^ffi{xi, k^,; A„ t,},=,i,..,n)- 

(27) 



IV. RESULTS AND DISCUSSION 



As we are interested in studying the effects of the meson cloud that notoriously manifest 
themselves at low values of Q^, the three free parameters of the model, /?, 7 in Eq. (|T3l) and 
the quark mass m, are fixed by fitting 8 experimental values of the proton and neutron form 
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Table II: Values of the form factors at = in the calculation of Ref. [Slij] (first column), from the 
bare nucleon contribution to the current matrix elements (second column) , from the contribution of 
the baryon (third column) and meson (fourth column) state in the BM component of the dressed 
nucleon. The column labeled TOT gives the total result in the meson-cloud model, while in the 
last column (labeled Exp.)are given the experimental values. 





Ref. [51] 






y- ^t{M'M)B 


TOT 


Exp. [62] 




2.78 


2.52 


0.18 


0.17 


2.87 


2.793 




-1.69 


-1.51 


-0.12 


-0.17 


-1.80 


-1.913 


gA 


1.24 


1.12 


0.075 


0.002 


1.20 


1.2670 



factors at low Q^. In the fit procedure we have used the Sachs form factors defined in terms 
of Dirac and Pauli form factors as 



AMI 



Gm{Q') = F,{Q')+F2{Q'). 



(28) 



The electric form factors are normalized as usual, i.e. 6*^(0) = 1, G'^(O) = 0, and the 
magnetic form factors at = are normalized to the nucleon magnetic moments, i.e. 
G^m{0) = /i^'". We have chosen to fit /i^, /i", the proton axial coupling constant = G'a(O), 
at g2 = 0.15 GeV^, and and G^/ at = 0.15 and 0.45 GeV^. 
A 5% uncertainty was allowed in the fitting procedure. The multidimensional integration 
required for the numerical computation was implemented in a parallel computation using 
the parallelized version of the VEGAS routine of Ref. 61]. 

The fitted values are 7 = 3.21, (3 = 0.489 GeV and m = 0.264 GeV. These values 
differ from the original set of parameters in Ref. 5l| , used for the calculation of the nucleon 
electromagnetic form factors in a three- valence-quark model of the nucleon. They were fitted 
only to the anomalous magnetic moments of the proton and neutron to obtain 7 = 3.5, 
p = 0.607 GeV and m = 0.263 GeV. ^ 

In TablellTlthe values of fi^, /x" and found in Ref. [51|] are compared with those obtained 



here and the experimental values 



formalism (see, e.g., Refs. 



37 



60 



As in Ref. 



63, 



64 



5l| and quite generally in the light-front 



65|). it is always difficult to reproduce all the 



three quantities simultaneously, so that some compromise has to be accepted. A significant 
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improvement is obtained for /i" by taking into account the meson cloud. One may appreciate 
that the two contributions in the baryon-meson fluctuation with the active particle being 
a meson {Y.b,b',m ^l'^^'^^^^) ^ baryon (Z^m.a/'.s ^-^'■*^'*^^'^) &dd up coherently in the right 
direction bringing the values of fi^, fi"' and qa closer to experiment. 




Figure 2: The four nucleoli electromagnetic form factors compared with the world data consid- 
ered in the analysis of Ref. [26] and the recent JLab data [19] using = + 
(5^/0.71GeV^)^ (open squares). Dotted (dashed) line for the contribution of the meson cloud 
(valence quarks). Solid line for the total result. 

The resulting electromagnetic form factors of both the proton and neutron are shown in 



Fig.[2]in comparison with the world data considered in the analysis of Ref. 26|] and the recent 



JLab data [19|. A rather good fit is obtained in the proton case in the whole range of available 
data, while in the neutron case the fit is less satisfactory. In any case, the contribution from 
the meson cloud is smooth and mainly significant for < 0.5 GeV^ with a maximum at 
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= 0. Therefore, in agreement with dispersion relation analyses 27|, |28|, |29| this model 
is unable to produce the bump/dip structure advocated in Ref. |26j around = 0.2 — 0.3 
GeV^. At higher Q"^ the explicit meson cloud contribution dies out, but indirectly affects the 



3J], which is equal to 0.91 in 



bare-nucleon contribution through the normalization factor Z 
our calculation. 

The results plotted in Fig. [2] are all obtained starting from an instant-form wavefunction 
"^^'^^^ in Eq. 0121) containing a totally symmetric S-wave part in the quark momenta. In 
order to improve the result for the neutron electric form factor it is known that the presence 
of a small admixture (1-2%) of mixed-symmetry 5"-wave components is most important 
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66, 



68l . l69l . l70l |. Following [69] we assume the mixed-symmetry S"-wave component to be 
represented by appropriate combination of mixed-symmetry spin-isospin wavefunctions with 
two radial wavefunctions of mixed symmetry of the form 



p,qj 



M4^vl;^'[^k{k.}.=i,3), 



^:''^kp,q) = Ar.^;^vl/A''W({k,W,3), (29) 



_^ qz pz _^ qz 

where \l/^'['^l({kj}j=i^3) is the symmetric 5*- wavefunction (JT3l) . Ms and A/'a are normalization 
factors, and p and q are the Jacobi coordinates 



P=-^^i^i + ^2), q= W-(ki-k2). 



(30) 



The consequences of including a small percentage of such mixed-symmetric contribution to 
the neutron electric form factor are illustrated in Fig. [31 Even a percentage as small as 1% 
is able to produce a quite good result compared to data. As anticipated in Ref. |69| the 
same calculation leaves the other nucleon form factors almost unaffected. 

The slope of the electric form factor at = determines the nucleon charge radius, i.e. 



-6 



dQ2 



(31) 



The corresponding values for proton and neutron obtained using an SU(6) symmetric or a 
mixed-symmetric instant-form wavefunction are reported in Table UTTt where also the partial 
contributions are indicated when one considers either the bare nucleon or the contribution of 
the cloud with an active baryon or meson. Rather good values of and r„ are obtained in 
the latter case. For the proton the charge radius is mostly due to the valence quarks, i.e. the 
bare proton. The meson cloud brings a contribution of about 5% which leads to a final value 
of 0.877 fm, in close agreement with the experimental value 0.8750 ± 0.0068 [62]. Including 
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Figure 3: The electric form factor of the neutron. Data points and sohd hne as in Fig. [2l Dashed 
(dotted) hne with 1% (2%) mixed-symmetry 5'-state in the bare neutron wavefunction. 

the mixed-symmetry S"-wave component the charge radius of the bare proton is 0.837 fm. 
The cloud adds a small contribution, that makes the total value to be slightly overestimated. 
In the neutron case the bare contribution is quite small, as expected. In contrast, both the 
cloud and the mixed symmetry are equally important. The contribution of the active meson 
in the cloud is substantial and with the right sign. Including also the mixed-symm etry 
S"-wave component becomes quite close to the experimental value —0.1161 ± 0.0022 62|. 

In Fig. m the electric to magnetic form factor ratio /i^G^/G^ is shown for the proton. 
The model follows the observed trend of data taken in polarized elastic electron scattering 
with a steepest fall-off at values of much larger than those involved in the fit to determine 
the model parameters. This is due to the combined effect of a slightly overestimated 
and a slightly underestimated at large values of (see Fig. [2]), where in practice only 
the valence quarks contribute. The result is only slightly modified when including a mixed- 
symmetry contribution in the nucleon wavefunction. 

The predicted axial form factor of the proton is shown in Fig. [51 The axial form factor 
has been normalized by its value at = 0, i.e. by the fitted value of the axial coupling 
constant qa- Also in this case the meson-cloud contribution is only significant at low values 
of Q^, although not sufficient to bring in complete agreement with experiment. However, 
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Figure 4: The electric to magnetic form factor ratio fi^G^/G^j for the proton. Data points as in 
Fig. [21 

the observed dipole form of the axial form factor, i.e. G aIQ"^) / G a{0) = 1/(1 + Q^/M|)^ 
with Ma = 1.069 GeV, is well reproduced. 

Neglecting relativistic corrections, in the Breit frame the radial distribution of the nucleon 
charge [pch) and magnetization (pm) are given by the Fourier- Bessel transform of the nucleon 
electromagnetic Sachs form factors, i.e. 

pTir) = - [ dQQ'jo{Qr)GnQ'), (32) 
f^^'^'pTir) = - I dQQ'jo{Qr)G%-{Q'). (33) 

71 J 

The corresponding results are shown in Fig. [6] together with the partial contributions from 
the bare nucleon and when the baryon or the meson is active in the cloud. In all cases the 
meson-cloud contribution is rather smooth and dies out beyond 2 fm. With the exception of 
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Table III: The different contributions (in fm^) to the proton and neutron mean square charge radii, 
Vp and respectively, from the bare nucleon and when the baryon or the meson is active in the 
cloud. The column labeled TOT is the total result. The lines labeled SU(6) (mixed symmetry) 
refer to the symmetry of the instant-form wavefunction of the bare nucleon. 







bare nucleon 


active baryon 


active meson 


TOT 


p 


SU(6) 


0.64 


0.065 


0.061 


0.77 


' p 


mixed symmetry 


0.70 


0.065 


0.061 


0.82 


' n 


SU(6) 


-0.0097 


0.0085 


-0.063 


-0.064 




mixed symmetry 


-0.058 


0.0085 


-0.063 


-0.112 



the neutron charge density, when the baryon is active its contribution is peaked at approxi- 
mately the same position as the bare-nucleon contribution, while the active meson is peaked 
at ~ 1 fm. Thus the meson cloud manifests itself as a slig ht extension of the radial distribu- 
tion up to ~ 2 fm as suggested by the analysis of Ref. 2a]. In the case of the neutron charge 
distribution the two components of the meson cloud behave differently. The active baryon, 
a proton or a A"*", gives a positive contribution of the same shape as the bare proton and is 
appropriately scaled by the corresponding vertex functions. The contribution of the active 
meson, a vr", is opposite and peaked at ~ 1.3 fm. The resulting charge distribution shows a 
positive core surrounded by a negative surface charge pushed outwards by the meson cloud 
and peaking at ~ 0.8 fm, in agreement with the analysis of Ref. 25|] and the expectation 



7l|. 



based on the picture of a hadron's periphery caused by the pion cloud 

The effect of including the mixed-symmetry S"-wave component in the neutron case can 
be appreciated from Fig. [71 The inner positive core is more pronounced and the negative 
surface charge is even more pushed outwards. Consequently, the (negative) mean square 
radius approaches the experimental value (see Table UTTl) . 

The direct relationship between Sachs form factors and the static charge and magnetiza- 
tion densities is lost when relativity is considered because there is a variation with of the 
Breit frame and electron scattering measures transitions matrix elements between nucleon 
states that have different momenta. Therefore one has to apply appropriate boosts that 
in a relativistic composite system such as the nucleon depend on the interaction among its 



18 




Figure 5: The axial form factor of the proton. Dotted (dashed) hne for the contribution of the meson 
cloud (valence quarks). Solid line for the total result. Dot-dashed line for the phenomenological 
dipole form. Data points are the world data considered in Ref. 7^ . 

constituents. The problem of finding a suitable prescription to relate Sachs form factors 



25| taking 



to the static charge and magnetization densities was recently addressed in Ref. 
into account the Lorentz contraction of the densities in the Breit frame relative to the rest 
frame. The consequences are that relativity tends to pull the density inward and to amplify 



oscillations at large radii. This has been confirmed in the model of Ref. [38|. The same 
effect should be expected also here. In any case in the present analysis the meson cloud is 
responsible for a long-range contribution to the nucleon charge and magnetization densities. 

The problem of unambiguously determining the charge density can be solved by looking at 
the charge density p{b) of partons in the transverse (impact parameter) plane with respect to 



the direction of the three-momentum transfer 30[ • This is possible because in the transverse 
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01230123 
r (fm) r (fm) 



Figure 6: The proton and neutron charge and magnetization densities. Dashed, dot-dashed and 
dotted Unes for contributions from the bare nucleon, the active baryon and the active meson in the 
cloud. SoUd hue for the total result. 



plane boosts are purely kinematica 



subgroup of the Poincare group [73 
of the Dirac form factor Fi : 



ji.e. in the light-front framework they form a Galilei 



7J]. Then p{b) is the two-dimensional Fourier transform 



1 f°° 



(34) 



where Jo is a cylindrical Bessel function. 

The corresponding charge densities for the proton and the neutron are plotted in Fig. [HI 



As in Ref. 



301 ] the densities are concentrated at low values of b with a positive peak for the 



proton and a negative peak for the neutron. 

These nucleon charge densities can be related to quark transverse distributions. Assuming 
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Figure 7: The neutron charge density with a mixed-symmetry S'-wave component in the neutron 
wavefunction. Dashed, dot-dashed and dotted hues for contributions from the bare nucleon, the 
active baryon and the active meson in the cloud. Sohd hne for the total result. 

that only up and down quarks are in the nucleon and invoking isospin symmetry, we have 

ff{b) = ^u{b) - ^d{b), (35) 
p"(6) = ~u{b) + ^d{b), (36) 

where u{b) is the transverse distribution for an up quark in the proton or a down quark in 
the neutron, and d{b) is the transverse distribution for a down quark in the proton or an up 
quark in the neutron. Both u{b) and d{b) are normalized to 1. They can be obtained using 

u{b) = fP{b) + ^p"(6), (37) 
d{b) = pf(6) + 2p'*(6). (38) 
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Figure 8: The proton and neutron charge density as a function of the impact parameter b. SoUd 
Unes for a permutationally symmetric momentum wavefunction, dashed hues with mixed-symmetry 
components included. 

The resulting distributions are shown in Fig. [9] in the two cases of a permutationally 
symmetric momentum wavefunction of the bare nucleon and of an included mixed-symmetric 
component. The central up quark density turns out to be larger than that of the down quark 
by about 40% in the symmetric case and by about 25% including the mixed symmetry. 



Quite similar results (about 30%) have been obtained in Ref. [30|| using phenomenological 
parametrizations of the Sachs form factors and deducing Fi in terms of G^; and Gm- 

The probability p'^{b, A, A) to find a quark with transverse position b and light-cone helicity 
A (= ±1) in the nucleon with longitudinal polarization A (= ±1) can be expressed as 
Fourier transform of the combination of the quark contributions to the Dirac and axial form 
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Figure 9: Transverse distributions of up (solid lines) and down (dashed lines) quarks in the proton 
as a function of the impact parameter b with a permutationally symmetric momentum wavefunction 
of the bare nucleon (left panel), and with a mixed-symmetric component (right panel). 



factors [731, i.e. 



AAG^(Q' 



2 ^ 



J dQQMQb) [Ff (g2) + AAG^(g2^ 
i [p^{b) + XAAq{b)] 



(39) 



where p'^{b) was already defined in Eq. flMl) and Aq{b) is the Fourier transform of (j^(Q^). 
Assuming a positive proton helicity (A = 1) the resulting probability is shown in the right 
panels of Fig. [101 The axial contributions Au{b) and Ad{b) for up and down quarks (left 
panels), respectively, have opposite sign. When suitably combined with the corresponding 
transverse distributions u{b) and d{b) given in Fig. [9] we see that the positive helicity up 
quarks in the proton are preferentially aligned with the proton helicity, while the opposite 
occurs for down quarks. This result is in total agreement with that shown in Fig. 7 of 
Ref. Its] where quite a different radial distribution of the axially symmetric spin density was 
presented for up and down quarks in the transverse plane. 
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Figure 10: Transverse distribution of up and down quarks in a longitudinally polarized proton as 
a function of the impact parameter b. Left panels: the axial contributions Au and Ad for up and 
down quarks, respectively. Right panels: total contribution for quarks polarized in the longitudinal 
direction, either parallel (solid lines) or antiparallel (dashed lines) to the proton helicity. 

V. CONCLUDING REMARKS 



The meson-cloud model, as revisited in Ref. [47] to study generalized parton distributions 
and including baryon-meson configurations with the baryon being a nucleon or a A and the 
meson being a pion as well as a vector meson such as the p or the u, has been used to describe 
the electroweak structure of the nucleon. Light-cone wavefunctions for the bare nucleon 



5l| and 



were constructed starting from the momentum wavefunction f[T3|) taken from Ref. 
depending on three parameters, the scale (3, the parameter 7 for the power-law behaviour, 
and the quark mass m. They are determined by fitting 8 experimental values of the proton 
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and neutron form factors at low Q^. No other free parameters enter the model calculations, 
since all other ingredients are fixed from the beginning on the basis of previous analysis. 

An overall good fit to the electromagnetic form factors is obtained, with the exception 
of the neutron electric form factor where it is essential to also include a mixed-symm etry 
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68, 



69|, 



7Q|. 



5"- wave momentum component, in agreement with earlier findings [66|, 
This component only slightly affect the other form factors. In any case the meson-cloud 
contribution is smooth and only significant below Q"^ = 0.5 GeV^. Therefore, as in analyses 



'27 



28 



based on dispersion relations 
the bump/dip structure discussed in Ref. 
proton axial form factor. 



29l |. also in this model no possibility exists to reproduce 
261. A similar smooth contribution arises in the 



Charge and magnetization densities are deduced as a function of both the radial distance 
from the nucleon center and the transverse distance (impact parameter) with respect to the 
direction of the three-momentum transfer. 

The meson cloud produces a slight extension of the radial distribution of the static charge 
and magnetization up to ~ 2 fm. It is confirmed that the neutron charge distribution shows 



30|. It is pushed outwards by 



a positive core surrounded by a negative surface charge 
the long-range meson cloud with opposite contributions from the active baryon and meson 
in the baryon-meson component of the nucleon wave function. 

As a function of the impact parameter a central negative charge is found for the neutron. 
A similar result has been obtained in Ref. 30| starting from a phenomenological fit of the 



electromagnetic form factors. This result can be explained invoking isospin symmetry and 
observing that the up quark transverse distribution in the proton is larger than the down 
quark one, a result consistent with deep inelastic scattering data. 

The up and down quark distributions associated with the Fourier transform of the axial 
form factor have opposite sign, with the consequence that the probability to find an up 
(down) quark with positive helicity is maximal when it is (anti)aligned with the proton 
helicity, in close agreement with the radial distribution of the axially symmetric spin density 



studied in Ref. 



75|. 



In conclusion, the meson-cloud model appears to give a satisfactory description of the 
electroweak properties of the nucleon with interesting information about its structure in the 
nonperturbative regime of QCD. 
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Appendix A: VERTEX FUNCTIONS 



In this Appendix we work out the case of the B B'V transitions in the light-front 
formalism, with the baryon states B being a nucleon or a A, and the vector mesons V 
corresponding to u, or p. The vertex functions for the coupling of baryons with pseudoscalar 
mesons are given in Appendix C of Ref . 47|] . The vertex functions for transitions to vector 



mesons can be found in several places (see, e.g., Refs. 



35 



and the longitudinal- 



momentum distributions corresponding to the integration over the transverse momentum of 

n n n 

the squared vertex functions are explicitly given in Refs. [73, [781, 179[. The quoted results 
are controversial, in the sense that they differ although the formalism is exactly the same. 



In particular, we agree with the resu 



Refs. 



78 



79], and differ from Ref. [77|, whi 



ts for the longitudinal-momentum distributions of 



conclusions drawn in the Appendix of Ref. 
Refs. 



35 



48 



br the vertex functions we agree with the 



79|], where the origin of the differences from 



54| is explained in details. However, Kumano et al. [79| do not give explicit 



analytical expressions for the vertex functions, and we find convenient to show here their 
derivation and in particular their dependence on the transverse momentum which enters in 
the convolution formulae for the form factors. 
The light-front vectors are defined as 



A'^ = {A-,A+,A^), 



(A.l) 



with 

We also use the notations Ar l = A^ ± iA^ and A = A_l 



(A.2) 
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The light-front nucleoli spinors Ux{p) are given by 



Ul/2{P) 



PR 
— m 

V ) 



U-l/2{p) 



( -p. ^ 

p^ + m 

PL 

ym — p^ J 



(A.3) 



The gamma matrices are defined as in Ref. 80 1. 

A similar expansion for the A field involves the Rarita-Schwinger spinors given by 



«3/2(P) 



iP) 



3/2 



e+i(p)wi/2(p), 

^eo(p)«-i/2(p) + ^e-i(p)Mi/2(p), 
e'ii(p)M_i/2(p), 



(A.4) 



where the polarization vectors are given by 

^'oiP) 



V2 v^' 



V P 

1 / 2 2 

1 / p]_ — m 



m \ p^ 



,P ,P± 



(A.5) 



The vertex function for the transition BV with the baryon B being one of the octet 
states is given by: 



9 u{Pb, ^'b)Y u{pn, A7v)e*,Av 



2M 



N 



(A.6) 



where F^bv is the vertex form factor, and the isospin factor is defined as 
{B\^*y ■f\N) = (-l)"''^^^=^(^^^ivl - ry\TBrB). 



(A.7) 



with Tb = I, and = ^6. 



27 



Table IV: Vertex functions for N N'V and particle helicities i — > A^, Ay in the prescription A. 



A'^ Av V{N,NV) 



,1 , 1 V2ki^ 



1-y ^ 2 



1,1 V2{MB^-yMN) / ^/2M4-(l-^/)2MJvMB] 



2 ^ y 



2Mjv 



'2 U y My^ 2Afjv v^(l-J/) 



1 _i / v^fej 

2 2Mjv v/l^(l-y) 



In the vertex function of Eq. (1A.6I) there is an off-shell dependence introduced by the 
derivative coupling, leading to a freedom in defining the vertex momentum. One can consider 



the following two possibilities 



54| 



(A) py = {Ev,pv), with Ev = ^m^+p^, 

(B) Py =pn -p'b = {En - E'b,pv). 



(A.8) 



In the following we will list the results for the vertex functions corresponding to both pre- 
scriptions, although the calculation of the nucleon form factors is performed using the off- 
shell condition (B), as suggested in Ref. [54 1. 

The results of the calculation for the particle helicities ^ A'^, Xy are given in Table HVl 
for prescription A and in Table IVl for prescription B. 

The corresponding results for helicity down of the nucleon are given by 

v-;ir'\y.kx) = (-i)''---i>f:'(!/,fcx), 

where k_L = {k^, —ky). 
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Table V: Vertex functions for N N'V and particle helicities i — > A^, Ay in the prescription B 



A'jv Ay V{N, NV) 



V2k 



Vv 



_9 I / 

1-y 2 



1,-1 ^/2{^IB-yMN) . f V2[kl-(MN+MB){l-y){yMM-MB)] 

fc2 ^(^l_y^2Mj^MB-yM^ 



g- 
+ 



2 " y Afvv^(l-J/) 

2Mn 2Mvy^ 



n ^ (Afa-Mjv) / kii{l+y)[k\-y(Ml+Ml+Ml)+Ml+y^M%] 
^ y Mv^ ^ 2Mm 2Mvy^(l-y) 



-lI _1 r, ^y^R I / V2kRMB 

^2 J- y Jy(i-y) ^ 2Mm Jy 



_l _i / 

2 ^ 2Afjv 



The vertex function for the transition BV with the baryon being one of the decuplet 
states is given by: 



f - ~ 

X Mp'b^ A') 757/. u{pn, A) [p^e^;; - 

(A.9) 



where the isospin factor is defined as in Eq. (jAlD, with = I and (|||T|||) = 2. 

The explicit results for particle helicities | — ^ A', Ay are given in Table IVTl for prescription 
A and in Table lylTl for prescription B. 

The corresponding results for helicity down of the nucleon are given by 
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Table VI: Vertex functions for N AV and particle helicities 5 — > A', Ay in the prescription A. 



A' Xv V{N, NV) 



,3 ,1 f 4 

^2 Mv y^(l- 



,3 n / MvVlkL 
-^2 ^ Mv Vvi^-y) 



I 3 _i / MNMsjl-yf-yMl 
"^2 ^ Mv. Wii-y) 



, 1 T 1 / kL{kl-2(l-y)M%) 
"•"2 VsMv MBy^/y(l-y) 



1 p, _ /2" / Mv[A:^-(l-i/)MB(MB-i/Mjv)] 
+ 2 U YsMv MBVy(l-2/) 



, 1 _i 1 / kR{-2MNMB{l-y)+yMl) 

"•"2 ^ a/3 MB^{l-y) 



_1 , J / 2MB(l-;/)fcj^-M|(l-j;)2+MjvM^i/3 

2 VsMv Msy^il-y) 



_1 n /2 / Mvfcfl[Mjvy-(l-j/)MB] 
2 ^ V 3 Ml, Ms ^^(1-2/) 



1 _i J / Mjvfc| 

2 ^ V3 M,^ Mb7j?(1-2/) 



3 -1 / MgkRjl-y) 
2 ^ Mv y^ 



-i 



-i -1 
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Table VII: Vertex functions for A'^ — > AV and particle helicities 5 A', Xy in the prescription B. 



A' Xv ViN, AV) 



^2 


+ 1 


/ 4 




^2 





f kLlk^,-yM^-y{l-y)M% + {l~y)M%] 
Mv V2Mvy^{l-y) 




^2 


-1 


f kl+{l-y)iMN+MB)(MB-yMN) 




Mv Vy{i-y) 




^2 


1 


1 / kL[kl-2(l-y)Ml] 

V3 Mv MByVy{l-y) 













V6 Mv MBMvy^{l-y) 




+ ^ 


-1 


1 / fcH[-fci+(l-2/)(2/M2,-2MivMs-M|)] 




y/zMv MB^{l-y) 




1 


1 


1 / fc2 [2Ms(l-j/)-2/2Mjv] + (l-y)[i/3W|,-MjvM2y2- 


M|(l-j/)] 


2 


V3 MBy^(l-y) 




1 





1 / fcfl[M;v2/-(l-2/)Ms][2/(l-2/)M2,-(l-y)M|+2/M2. 




2 


a/6 A^V' MBMvy^/y{l-y) 




1 
2 


-1 


1 / Mjvfel 
V3 MB^{l-y) 




3 
2 


1 


f MBkR{l-y) 
Mv yx/y 




3 
2 










3 
2 


-1 
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